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Employing an idea of Brekhovskikh, an expression for the reflection coefficient of a 
continuously stratified ionized medium is derived. The result is in the form of a series 
whose first term is a Fresnel-type coefficient and succeeding terms account for the finite 
thickness of the transition layer. This result is then fitted into previously developed theory 
for propagation between a spherical earth and a concentric ionosphere. 

1. Introduction 

In the mode theory of VLF radio propagation it is often assumed that tlie lower edge of 
the ionosphere may be represented by a sharply bounded and homogeneous ionized medium 
[Wait, I960]. This assumption could be criticized as being unrealistic as it is known that the 
boundary between tlie non-ionized lower atmosphere and the ionosphere is ''graduaL" Un- 
fortunately, the profiles of electron density and collision frequency, in the height range of 
importance for reflection of VLF waves at oblique incidence, are not well known. At these 
low heights (i.e., around 70 km), the rocket measurements are not reliable and vertical inci- 
dence ionosondes give even less information. Some important qualitative data have come from 
the use of weak backscattered echoes at medium frequencies [Ratcliffe, I960]. Among the 
findings of such work is the occurrence of very sharp gradients in electron density in tlie height 
regions of 60 to 80 km. 

In view of the preceding remarks, the sharply bounded model is not completely uujiistified. 
Furthermore, it agrees with much of the experimental data at VLF if attention is restricted to 
highly oblique incidence. Nevertheless, it is desirable if more attention is paid to the possible 
influence of the '^gradualness'^ of the ionosphere boundary. By making use of a digital com- 
puter some excellent progress has been made in this direction [Ferraro and Gibbons, 1959; 
Barron, 1959; Johler, 1962; Budden, 1961]. In particular, Johler [1961] has made a careful 
comparison between full wave solutions for certain gradual and sharpl}^ bounded models. His 
results indicate that the sharp boundary is often quite adequate. 

It is the purpose of this paper to outline a relatively simple approach to the problem and 
to fit the results into the mode theory. 

2. Formulation 

To start with we shall assume an inhomogeneously stratified medium. With respect to a 
Cartesian coordinate system (x,y,z), the refractive index N{z) is assumed to be only a function 
of 0. As z tends to positive and negative infinity the refractive index is assumed to approach 
constant values. That is 

and 

At 2=— 00^ a plane wave is incident at an angle ^o w^ith respect to the positive z direction. 



1 The research reported in this paper has been sponsored by the Advanced Research Projects Agency, Washington 25, D.C. 
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To simplify the discussion the earth's magnetic field is neglected. This is justified, at 
least in part, since the collision frequency is somewhat greater than the gyro frequency of 
electrons in the lowest ionosphere. 

Thus A^(^)=ri-i^T 

where co is the angular f requeue}^, and a;^(^)=— -V (1) 

where coo and jjl are the (angular) plasma and collision frequencies, respectively [Ratcliffe, 1959]. 
These latter quantities are also functions of the coordinate z. It is convenient to introduce a 
dimensionless factor Liz) in the following manner. 



^w=(-ife)' 



where 1 ^i^M^^i^ (2) 

L{Z) CO (v + lco)aj 

Then 

and 

In most applications at VLF, z^>>co, so that L{z) can be regarded as a real parameter. 

The problem is to calculate the reflected wave at z=~ oo. If the electric vector of the 
incident wave was perpendicular to the ^ axis, the resulting differential equations can be solved 
in terms of known functions provided certain profiles of N{z) are used. The elegant analysis 
of P. S. Epstein [1930] for normal incidence can be carried over directly to this case (i.e., for 
horizontal polarization). Other closed-form solutions for the case of ^'-parallel polarization 
are also available in the literature [Brekhovskikh, I960]. Unfortimately, the differential 
equations for i7-parallel polarization (i.e., for vertical polarization) are somewhat more com- 
plicated and closed-form solutions for special transition profiles do not seem to be readilj^ 
available. Thus a perturbation procedure is adopted following the suggestion of Brekhovskikh 
[I960]. 

3. Differential Equation for the Reflection Coefficient 

Without any loss of generality, the H vector is taken to have only a y component. Thus, 
for a time factor, e'^^^ MaxwelPs equations are given by 

-^ — ^ = -^Mcoi?„ iei^E,=^ ^^, and teco£;,=-^ (3) 

where fi and e are the magnetic permeability and permittivity of the medium. For sake of 
generality, both jj, and e can be regarded as functions of z. Of course, in applications to the 
ionosphere ju can be replaced by its free space value juq. The retention of a variable /x in the 
theory permits one to readily adapt the results to arbitrary polarization. Furthermore, the 
analog,y in acoustics is readily brought out. 

Now the sum of the incident and reflected waves is defined in the following manner. 
The sole magnetic field component is written 

Hy=[A{z)+B{z)] exp [-i{eiiy^Sx\ (4) 

while the two electric field components are written 

E,= [A{z)-B{z)] aM/€)*exp [-i{eixfSx] (5) 

54 



and 

E,= -[A{z)+B{z)]S{^ler^ exp [-^(6M)^&) (6) 

in terms of S and C^ wliich can be regarded as sine and cosine of angle 6 which is itself a function 

of z. At z=— CO the incident wave is taken to have an angle of incidence ^o. 

In tliese equations, /:==(eoiLto)'<^ where cq and mo axe the constants of free space, while A{z) and 

B{z) are not yet defined. Since these field components are to satisfy Maxwell's equations, it is 

required that 

(6M)^S'=constant 

and since S^sin 6q, e = eo, m=Moj at ^ = — oo^ it follows that 

{e^xy^S={e,ixo)Umd, (7) 

which is just a statement of Snell's law. Furthermore, 

Again, as a consequence of Maxwell's equations, A{z) and B{z) must satisfy 

^+i5A+V{A-B) = Q (8) 

^-i5B+r(B-^) = (9) 

dz 

where 
and 



r= 



2d dz \€/ 



Equations (8) and (9) are easily combined ijito a single equation for the ratio BjA = R{z). Thus 

^=2idR + T{l-R') (10) 

where R(z) is, by definition, a reflection coefRciejit. Results, more or less equivalejit to (10) 
have been given by Budden [1961] and Brekhovskikh [I960]. 

4. Iterative Solution 

To obtahi a solution, R{z) is written in terms of a new function v{z). Thus 

g{z)v(z)+gi 

Now, since ^^^ R{z) = 0, it follows that ^^^^ v{z) = l. The differential equation for v{z) 'is 
obtained by substituthig (11) into (10). This can be written in the relativel}^ simple form 

-i-T^('-i;-) <-' 

where A*=(€omo)^ ^ and A^, K, and g are functions of z. Usijig a method of successive approxi- 
mations (see for example, Z. Kopal [1955]), the solution can be expressed as an ascending 
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series in powers of k. For example, the zeroth approximation is to replace the right-hand 
side of (12) by zero, thus ?; is a constant which must be unity to satisfy the limiting condition 
at 2^ 00 . The first approximation is obtained by replacing v'^ on the right-hand side by unity. 
Thus 

v=l + ikg,y^ ^{l-^j^z (13) 

where the limits of the integration are chosen so that v satisfies the limiting condition at ^== oo . 
The second approximation is then obtained by substituting the latter result for v into the 
right-hand side of (12). In general, the nth approximation, Vn, can be found from the (n-l)th 
approximation, Vn-\j by using 

Vr,=l + ikg, J^" ^ (l-|^ vl.^dz (14) 

5. Some Simple Extensions of the Solution 

While these results have been developed with specific reference to an incident wave with 
the electric vector in the plane of incidence, the results are also applicable to the other polari- 
zation. If the magnetic vector of the incident wave is in the plane of incidence (i.e., horizontal 
polarization), the results are still valid if the following transformations are made. 

Hy-^Ey, E^-^—H^, E^-^—H^, fx->e, and e-^^. 

Thus the formula for the reflection coefliicient R, given by (11), is still valid if K is replaced 
by €/€o and N is not changed. 

There is also a well-defined acoustic analogy to the problem being discussed. In this case 

A^(^)=-T^ where c(z) is the velocity of sound and Cq is the limiting value of c at 2=— oo . Thus 

the velocity is varying from a constant value Cq to a differing constant value Ci at z=-{- oo. 
Also, K(2)=p(2)Ipo where p(2) is the density and po is its limiting value at 2= oo. The com- 
ponent Hy is then analogous to the acoustic pressure and Ex and Ez are analogous to the x and 2 
components, respectively, of the particle velocity. 

6. Discussion of the Form of the Solution 

It is interesting to note that for the zeroth approximation (corresponding to v=l), the 
reflection coefficient may be written 

R= T^ where Co=cos do and ^So^sin Oq. (15) 

This is the Fresnel reflection coefficient for the reflection of a plane wave at oblique incidence 
from a sharply boimded and homogeneous medium. Thus the higher terms in the ascending k 
series account for the ^'gradualness'^ of the boundary. 

In the general case, the reflection coefficient may be written 

^^^"^ Co+{Ki/Ni)im-si)yvizo) ^'''^ 

where Zo is some convenient level, below which K(z) and N{z) may be regarded as unity. Thus 
the total field in the region z<iZ(, can be written 

i7^=i7oe-*'^«o[e-'*<^-VCo4-_B(2„)e+«(-'+^o'Co] (17) 
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where Hq is the value of the incident wave at the fictitious interface, 2=Zq. In the second 
approximation, neglecting terms in F and higher, 

Succeeding terms quickly become more complicated. It can be seen that the integrands con- 
tain the factor l — ig/giY in each of these terms. The presence of this factor permits one to 
replace the upper limit of each of these integrals by Zi, where 2i is the level above which N{z) 
and K{2) may be replaced by Ni and Ki. Thus the transition region may be defined as the 
interval ^o<C^<C^i- On this basis, it is apparent that the nth term in the scries for v(3o) is the 
order of [k(3i — 2Q)]^. Consequently, the series converges rapidly when the electrical thickness 
of tlie transition layer is small. 

As an interesting check, the transition is replaced by a homogeneous slab. Tluis, 

N(z)=N and K{z)=K when 2o<^<^i. 
In this case, the reflection coefficient is given exactly by Wait [1958]. 



provided 



«^(^o) = 



j._ Cov(2o)-(Ni-SinKJN-i) . 

^'^A K ^^^^ [ik(zi—ZonN'-S'o] 

1_| i^ ^^ tanh [%k[Zi — Zo)yN'^—SQ] 



If |/:(^i— ^o)ViV^— /S'o|«l, it is seen that to a first order, 



>. 



»(2o)^l+ j^2 



l-^Ar 4-^^ \ik{z,-z,)^ (21) 



which is consistent with the first two terms of the series given by (13). 

For some applications to VLF propagation, it is desirable to express the reflection coeffi- 
cient in the following form 



where 






Then, if ^Cq is regarded as a small parameter, the following expansion results 

B=-e-'^'^+\ {^C,Y-\ {^C,Y+ . . . (23a) 

Thusif liSaP^Cl, 

Xl^_^-2^co (23b) 

which is a convenient form when the incidence is highly oblique and the frequency is not too 
low. On the other hand, if ^Cq is regarded as a large parameter, it is convenient to use the 
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expansion 

i?=exp {-£)-l (^3+1 (^.+ . . . (24) 

which can be approxmiated by the first term if \pCo\^^l. 

7. Application to Terrestrial VLF Propagation 

The results developed in the above form may be readily fitted into VLF propagation 
theory which has been developed recently in a very general form [Wait, 1961]. The source 
of the field is a vertical electric dipole located on the surface of a smooth spherical earth of 
radius a, conductivity a and dielectric constant e. Spherical coordinates are chosen with the 
dipole located at r=a and 6=0. A concentric reflecting shell is now located at r=a -\-h. The 
electrical properties of the interfaces at r=a and r=a-\-h are defined in terms of two dim en- 
sionless factors, q and qt, respectively. These are defined explicitly below\ Then, for har- 
monic time dependence, the electric field is written, apart from a constant factor, in the form 
[Wait, 1961] 

- ikad 



where 



Er= ,; . ^.. V (25) 



v= 2 y, (26) 

j= 0,1,2, . . . 



The Vj can be identified as sky waves or ^'hops'\ as they have a precise ray interpretation in 
the domain where geometrical optics is valid. Actually, in terms of the parameters of the 
problem [Wait, 1961], 

iov j even, while 

for j odd. In the above, x={kaj2)^^^ 6 and y={2lkaY'^ k{r—a). 

The contours of the integrations are to enclose all singularities of the integrands. The symbol 
t has the following meaning 

{-ty=(~)''c (28) 

where C can be geometrically interpreted as the cosine of the (complex) angle of incidence on 
the earth^s surface at r=a. Furthermore, 

(2/o-0*=(y)'"c" (29) 

where ^o=f — j kh and 6" can be geometrically interpreted as the cosine of the (complex) 

angle of incidence at the surface, r=a^h, which is regarded as the lower edge of the ionosphere. 
Other factors, found in (27), are 



and 



^(o=-[" ^;;;-^°;+g-^';;-y°n m 

U^i(t)-qwi(t)j 



where Wi (t) and W2{t) are Airy functions, defined elsewhere [Wait, 1961]. The factor q occurring 
in the above is defined by 
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since cr»eco in most cases of interest, this can be written 

/7.^,\l/3 



where 



G^-=— • 



To relate the factor g^ to the properties of the ionosphere reflecting layer, a somewhat 
heuristic approach is used. First, it is assumed that yo—t has a real part somewhat greater 
than unity; thus, Wi(t—yo), etc., may be replaced by the first term of their asymptotic expan- 
sions. Then 

^(„,g;^exp[-.:|-i|(|)<C').] (.34, 

/ 2\^/^ . 
where Ai=^lj—] iqt. As indicated in a previous analysis [Wait, 1961], (C — Ai)/(f ' + A^) is a 

reflection coefficient referred to the bottom of the layer r=a-\-h. Therefore, Aj may be identi- 
fied with the factor 1//3; this is certainly evident on comparing (22a) and (34). Implicit in the 
above identification is that the local reflection process at the ionosphere does not depend on 
eartli curvature. It slioukl be emphasized that, in general, ([t and (3 are functions of (^^ and 
thus, in turn, they are functions of /, the variable of integration. However, for propagation to 
great distances, q^ and may l)o rej)hi('ed by a (complex) constant appropriate for grazing in- 
cidence (i.e., So^l). 

On utihzing (2), (21), and (22b), it is seen that to a first order in kj 



or 



1 

% F 



^"1 7-2 



(36) 



Then, iqi=(ka/2y^^ (1/iS). In the above, it is understood that L is a function of 2 in the tran- 
sition region extending from 2q to z. If k(zi— Zq) is not small compared with unity, it is necessary 
to obtain higher order terms. 

As has been demonstrated before [Wait, 1961], the waveguide modes are found from the 
equation 

A(0S(0 = 6-^'2"« (37) 

where n is an integer. Some of the methods developed for solving this equation in the sharply 
bounded limit (i. e., Zi-^Zq) are also applicable to the '^graduaP' model. 

8. Application to ELF Propagation 

At ELF (extremely low frequencies), the influence of the finite transition region can be 
quite sun ply indicated. In this case, the niodal equation can be written 

B,{C)Ri{C)e-'''''^=e-'^''^ (38) 

where 

i?,(C)^exp[-(2/«7)Vi^] (39) 

and 

iZ,(C)^expl-(2/6')(l/^)]. - (40) 
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These are valid subject to the restrictions that iC^|>> (-)> vF")^' I'^l^l^ ^^^^ V^^- Equation 
(38) may be written 



.A 

■n-\- 

This is solved as a quadratic to yield 



khC=7rn+i -F^ where A=V^6''+1//?. 



The positive sign before the radical is chosen since it reduces to Cn= (irn/kh) when A approaches 
zero. The corresponding values of t are then found from 

Individual modes propagate in the horizontal direction according to (sin 6)~-' exp (—ixtn) 
exp {—ikad) which is equivalent to 

[sin (d/a)]~^ exp (—ikdSn), with d=ad, 

where Sn=(l — C'i)^. Since |A| fcA<<l, the radical in (41) may be expanded to yield 

Usually the zero-order mode is of greatest interest since the higher modes are beyond ^'cut-off.'' 
Furthermore, in most cases |A|<^<^M; therefore, the propagation is adequately described by 
the zero-order mode which may be expressed in the form 

The attenuation in nepers per unit distance is thus given by 

-ArImS'„^^[Re^+V^] (44) 



and 



where 






When b tends to zero, corresponding to a sharp boundary, tlie integral term vanishes and the 
attenuation factor is given by 

-kImSo^-^[^i+M 

which is well-known. When the dependence of i(^) with height is specified it is a relatively 
simple matter to calculate the correction resulting from the ^^gradualness" of the boundary. 
In general it can be seen that, for the ELF range, the contribution from the second integral 
is always positive since i>>Li over the transition region. Therefore, the presence of the 
transition layer is to increase the attenuation relative to that of a sharply bounded model. 

Specific numerical results for mode characteristics based on this approach will be given 
in a later paper. In particular, it is hoped to make use of the extensive numerical results of 
Johler [1961, 1962] for the plane wave reflection coefficients. 
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